Incorporating patient-specific hip orientation from weightbearing computed tomography affects discrete element analysis-computed regional joint contact mechanics in individuals treated with periacetabular osteotomy for hip dysplasia.

Computational models of the hip often omit patient-specific functional orientation when placing imaging-derived bony geometry into anatomic landmark-based coordinate systems for application of joint loading schemes. The purpose of this study was to determine if this omission meaningfully alters computed contact mechanics. Discrete element analysis models were created from non-weightbearing (NWB) clinical CT scans of 10 hip dysplasia patients (11 hips) and oriented in the International Society of Biomechanics (ISB) coordinate system (NWB-ISB). Three additional models were generated for each hip by adding patient-specific stance information obtained via weightbearing CT (WBCT) to each ISB-oriented model: (1) patient-specific sagittal tilt added (WBCT-sagittal), (2) coronal and axial rotation from optical motion capture added to (1; WBCT-combo), and (3) WBCT-derived axial, sagittal, and coronal rotation added to (1; WBCT-original). Identical gait cycle loading was applied to all models for a given hip, and computed contact stress and contact area were compared between model initialization techniques. Addition of sagittal tilt did not significantly change whole-joint peak (p = 0.922) or mean (p = 0.871) contact stress or contact area (p = 0.638). Inclusion of motion-captured coronal and axial rotation (WBCT-combo) decreased peak contact stress (p = 0.014) and slightly increased average contact area (p = 0.071) from WBCT-sagittal models. Including all WBCT-derived rotations (WBCT-original) further reduced computed peak contact stress (p = 0.001) and significantly increased contact area (p = 0.001). Variably significant differences (p = 0.001-1.0) in patient-specific acetabular subregion mechanics indicate the importance of functional orientation incorporation for modeling applications in which local contact mechanics are of interest.

The Effects of Residual Femoral Deformity on Computed Contact Mechanics in Patients Treated With In Situ Fixation for Slipped Capital Femoral Epiphysis.

OBJECTIVE: In situ fixation for treatment of slipped capital femoral epiphysis (SCFE) can stabilize the epiphysis and prevent further joint deformation but often leaves residual deformity that may adversely affect intra-articular contact mechanics. The purpose of this study was to investigate the relationship between residual deformity and contact mechanics in the post-SCFE hip. METHODS: Patient-specific hip models were created for 19 patients with SCFE treated with in situ fixation. For each model, discrete element analysis was used to compute cumulative acetabular and femoral contact stress exposure during a walking gait cycle. Slip severity was evaluated for each patient using the two-dimensional Southwick angle and a novel three-dimensional (3D) assessment of multiplanar femoral deformity (3D slip angle). RESULTS: Of the SCFE cases, 2/7 mild (Southwick angle ≤30 degrees) had peak cumulative femoral exposures equivalent to that of severe (Southwick angle ≥60 degrees) cases. Severe SCFE cases had higher peak ( P = 0.015) and mean ( P = 0.028) femoral contact stress exposure and lower cumulative femoral contact area ( P = 0.003) than mild (Southwick angle ≤30 degrees) SCFE cases. Mean femoral contact stress exposure was also higher in severe SCFE cases than in moderate SCFE cases ( P = 0.027). Acetabular and femoral contact mechanics metrics typically demonstrated stronger correlations with 3D slip angle than two-dimensional Southwick angle. CONCLUSIONS: Increased slip severity adversely impacts intra-articular femoral contact mechanics. Contact mechanics metrics demonstrate higher correlations with 3D slip angle, indicating that this novel measurement may better describe global deformity and its relationship to intra-articular mechanics; however, the modest strength of these correlations may also imply that global impingement-generating deformity is not the primary factor driving contact mechanics in the post-SCFE hip. CLINICAL RELEVANCE: Greater slip severity adversely impacts contact mechanics in the post-SCFE hip. However, focal regions of high contact stress were seen even in mild SCFE deformities, suggesting some type of deformity correction should be considered even for mild slips to alleviate secondary impingement, address focal incongruities, and reduce osteoarthritis development/progression.

Providing a computationally derived, mechanically optimised target correction during preoperative planning can improve joint contact mechanics of hip dysplasia treated with periacetabular osteotomy.

AIM: Preoperative identification of acetabular corrections that optimally improve joint stability and reduce elevated contact stresses could further reduce osteoarthritis progression in patients with hip dysplasia who are treated with periacetabular osteotomy (PAO). The purpose of this study was to investigate how providing patient-specific, mechanically optimal acetabular reorientations to the surgeon during preoperative planning affected the surgically achieved correction. METHODS: Preoperative CT scans were used to create patient-specific hip models for 6 patients scheduled for PAO. A simulated acetabular fragment was extracted from the preoperative pelvis model and computationally rotated to simulate candidate acetabular reorientations. For each candidate, discrete element analysis was used to compute contact stresses during walking, which were summed over the gait cycle and scaled by patient age to obtain chronic contact stress-time exposure. The ideal patient-specific reorientation was identified using a cost function that balances minimising chronic stress exposures and achieving surgically acceptable acetabular coverage angles. The optimal reorientation angles and associated contact mechanics were provided to the surgeon preoperatively. After PAO was performed, a model of the surgically achieved correction was created from a postoperative CT scan. Radiographic coverage and contact mechanics were compared between preoperative, optimal, and surgically achieved orientations. RESULTS: While surgically achieved reorientations were not significantly different from optimal reorientations in radiographically measured lateral (p = 0.094) or anterior (p = 0.063) coverage, surgically achieved reorientations had significantly (p = 0.031) reduced total contact area compared to optimal reorientations. The difference in lateral coverage and peak chronic exposure between surgically achieved and optimal reorientations decreased with increasing surgeon experience using the models (R² = 0.758, R2 = 0.630, respectively). CONCLUSIONS: Providing hip surgeons with a patient-specific, computationally optimal reorientation during preoperative planning may improve contact mechanics after PAO, which may help reduce osteoarthritis progression in patients with hip dysplasia.

Radiographically successful periacetabular osteotomy does not achieve optimal contact mechanics in dysplastic hips.

BACKGROUND: Optimal correction of hip dysplasia via periacetabular osteotomy may reduce osteoarthritis development by reducing damaging contact stress. The objective of this study was to computationally determine if patient-specific acetabular corrections that optimize contact mechanics can improve upon contact mechanics resulting from clinically successful, surgically achieved corrections. METHODS: Preoperative and postoperative hip models were retrospectively created from CT scans of 20 dysplasia patients treated with periacetabular osteotomy. A digitally extracted acetabular fragment was computationally rotated in 2-degree increments around anteroposterior and oblique axes to simulate candidate acetabular reorientations. From discrete element analysis of each patient's set of candidate reorientation models, a mechanically optimal reorientation that minimized chronic contact stress exposure and a clinically optimal reorientation that balanced improving mechanics with surgically acceptable acetabular coverage angles was selected. Radiographic coverage, contact area, peak/mean contact stress, and peak/mean chronic exposure were compared between mechanically optimal, clinically optimal, and surgically achieved orientations. FINDINGS: Compared to actual surgical corrections, computationally derived mechanically/clinically optimal reorientations had a median[IQR] 13[4-16]/8[3-12] degrees and 16[6-26]/10[3-16] degrees more lateral and anterior coverage, respectively. Mechanically/clinically optimal reorientations had 212[143-353]/217[111-280] mm2 more contact area and 8.2[5.8-11.1]/6.4[4.5-9.3] MPa lower peak contact stresses than surgical corrections. Chronic metrics demonstrated similar findings (p ≤ 0.003 for all comparisons). INTERPRETATION: Computationally selected orientations achieved a greater mechanical improvement than surgically achieved corrections; however, many predicted corrections would be considered acetabular over-coverage. Identifying patient-specific corrections that balance optimizing mechanics with clinical constraints will be necessary to reduce the risk of osteoarthritis progression after periacetabular osteotomy.

Joint contact stress improves in dysplastic hips after periacetabular osteotomy but remains higher than in normal hips.

AIM: The purpose of this study was to use computational modeling to determine if surgical correction of hip dysplasia restores hip contact mechanics to those of asymptomatic, radiographically normal hips. METHODS: Discrete element analysis (DEA) was used to compute joint contact stresses during the stance phase of normal walking gait for 10 individuals with radiographically normal, asymptomatic hips and 10 age- and weight-matched patients with acetabular dysplasia who underwent periacetabular osteotomy (PAO). RESULTS: Mean and peak contact stresses were higher (p < 0.001 and p = 0.036, respectively) in the dysplastic hips than in the matched normal hips. PAO normalised standard radiographic measurements and medialised the location of computed contact stress within the joint. Mean contact stress computed in dysplastic hips throughout the stance phase of gait (median 5.5 MPa, [IQR 3.9-6.1 MPa]) did not significantly decrease after PAO (3.7 MPa, [IQR 3.2-4.8]; p = 0.109) and remained significantly (p < 0.001) elevated compared to radiographically normal hips (2.4 MPa, [IQR 2.2-2.8 MPa]). Peak contact stress demonstrated a similar trend. Joint contact area during the stance phase of gait in the dysplastic hips increased significantly (p = 0.036) after PAO from 395 mm2 (IQR 378-496 mm2) to 595 mm2 (IQR 474-660 mm2), but remained significantly smaller (p = 0.001) than that for radiographically normal hips (median 1120 mm2, IQR 853-1444 mm2). CONCLUSIONS: While contact mechanics in dysplastic hips more closely resembled those of normal hips after PAO, the elevated contact stresses and smaller contact areas remaining after PAO indicate ongoing mechanical abnormalities should be expected even after radiographically successful surgical correction.

Effect of modeling femoral version and head-neck offset correction on computed contact mechanics in dysplastic hips treated with periacetabular osteotomy.

While correction of dysplastic acetabular deformity has been a focus of both clinical treatment and research, concurrent femoral deformities have only more recently received serious attention. The purpose of this study was to determine how including abnormalities in femoral head-neck offset and femoral version alter computationally derived contact stresses in patients with combined dysplasia and femoroacetabular impingement (FAI). Hip models with patient-specific bony anatomy were created from preoperative and postoperative CT scans of 20 hips treated with periacetabular osteotomy and femoral osteochondroplasty. To simulate performing only a PAO, a third model was created combining each patient's postoperative pelvis and preoperative femur geometry. These three models were initialized with the femur in two starting orientations: (1) standardized template orientation, and (2) using patient-specific anatomic landmarks. Hip contact stresses were computed in all 6 model sets during an average dysplastic gait cycle, an average FAI gait cycle, and an average stand-to-sit activity using discrete element analysis. No significant differences in peak contact stress (p = 0.190 to 1), mean contact stress (p = 0.273 to 1), or mean contact area (p = 0.050 to 1) were identified during any loading activity based on femoral alignment technique or inclusion of femoral osteochondroplasty. These findings suggest that presence of abnormal femoral version and/or head-neck offset deformities are not themselves predominant factors in intra-articular contact mechanics during gait and stand-to-sit activities. Inclusion of modified movement patterns caused by these femoral deformities may be necessary for models to adequately capture the mechanical effects of these clinically recognized risk factors for negative outcomes.

Chronically elevated contact stress exposure correlates with intra-articular cartilage degeneration in patients with concurrent acetabular dysplasia and femoroacetabular impingement.

Hip dysplasia is known to lead to premature osteoarthritis. Computational models of joint mechanics have documented elevated contact stresses in dysplastic hips, but elevated stress has not been directly associated with regional cartilage degeneration. The purpose of this study was to determine if a relationship exists between elevated contact stress and intra-articular cartilage damage in patients with symptomatic dysplasia and femoroacetabular impingement. Discrete element analysis was used to compute hip contact stresses during the stance phase of walking gait for 15 patients diagnosed with acetabular dysplasia and femoral head-neck offset deformity. Contact stresses were summed over the duration of the walking gait cycle and then scaled by patient age to obtain a measure of chronic cartilage contact stress exposure. Linear regression analysis was used to evaluate the relationship between contact stress exposure and cartilage damage in each of six acetabular subregions that had been evaluated arthroscopically for cartilage damage at the time of surgical intervention. A significant correlation (R2 = 0.423, p < 0.001) was identified between chondromalacia grade and chronic stress-time exposure above both a 1 MPa damage threshold and a 2 MPa-years accumulated damage threshold. Furthermore, an over-exposure threshold of 15% regional contact area exceeding the 1 and 2 MPa-years threshold values resulted in correct identification of cartilage damage in 83.3% (55/66) of the acetabular subregions loaded during gait. These results suggest corrective surgery to alleviate impingement and reduce chronic contact stress exposures below these damage-inducing thresholds could mitigate further cartilage damage in patients with hip dysplasia.

Isolated changes in femoral version do not alter intra-articular contact mechanics in cadaveric hips.

Abnormal femoral version is a deformity in the angle between the femoral neck and the transcondylar axis of the knee. Both femoral anteversion and retroversion alter passive and active rotation of the hip and are associated with intra-articular or extra-articular impingement. However, little is known about the effect of abnormal femoral version on intra-articular hip contact stresses. To quantify the effect of femoral version on hip contact stress, five cadaveric pelvis specimens were mechanically tested with a hip-specific Tekscan sensor inserted in the joint space. Specimens were oriented in a heel-strike position and loaded with 1000 N of compressive force. Pressure measurements were recorded by the Tekscan sensor with the femur oriented in 0°, 15°, and 30° of version. At the completion of testing, specimens were locked into place at 0° and post-test CT scans were obtained to register the pressure sensor measurements to the joint anatomy. There were minor changes in contact area (<7%) and translation of the peak contact stress location (8.8 ± 7.6 mm). There was no significant change in peak contact stress (p = 0.901) in either the retroverted (0°) or anteverted (30°) conditions relative to normal version (15°) under identical gait-related loading conditions. While abnormalities in patient gait and resultant joint loading caused by femoral version abnormalities may contribute to hip pain, the present findings would suggest that future joint degeneration in hips with version abnormalities are not simply the result of abnormal contact stress induced by joint incongruity due to femoral version abnormalities.

Patient Age and Hip Morphology Alter Joint Mechanics in Computational Models of Patients With Hip Dysplasia.

BACKGROUND: Older patients (> 30 years) undergoing periacetabular osteotomy (PAO) to delay THA often have inferior patient-reported outcomes than younger adult patients (< 30 years). It is unclear how patient age affects hip morphology, mechanics, or patient-reported outcome scores. QUESTIONS/PURPOSES: (1) Is increased patient age associated with computationally derived elevations in joint contact stresses? (2) Does hip shape affect computationally derived joint contact stresses? (3) Do computationally derived joint contact stresses correlate with visual analog scale (VAS) pain scores evaluated at rest in the clinic at a minimum of 1 year after surgery? METHODS: A minimum of 1 year of clinical followup was required for inclusion. The first 15 patients younger than 30 years of age, and the first 15 patients older than 30 years of age, who underwent PAO for treatment of classic dysplasia (lateral center-edge angle < 25°) who met the minimum followup were selected from a historical database of patients treated by a single surgeon between April 2003 and April 2010. The older cohort consisted of 14 females and one male with a median age of 41 years (range, 31-54 years). The younger cohort consisted of 10 females and five males with a median age of 19 years (range, 12-29 years). Median followup for the older than 30 years versus younger than 30 years cohort was 19 months (range, 12-37 months) versus 24 months (range, 13-38 months). Pre- and postoperative hip models were created from CT scans for discrete element analysis (DEA) contact stress computations. DEA treats contacting articular surfaces as rigid bodies (bones) separated by a bed of compressive springs (cartilage), the deformation of which governs computation of joint contact stresses. This technique greatly simplifies computational complexity compared with other modeling techniques, which permits patient-specific modeling of larger cohorts. Articular surface shape was assessed by total root mean square deviation of each patient's acetabular and femoral cartilage geometry from sphericity. Preoperative and postoperative VAS pain scores evaluated at rest in the clinic were correlated with computed contact stresses. RESULTS: Patients older than 30 years had higher predicted median peak contact stress preoperatively (13 MPa [range, 9-23 MPa; 95% confidence interval {CI}, 11-15 MPa] versus 7 MPa [range, 6-14 MPa; 95% CI, 6-8 MPa], p < 0.001) but not postoperatively (10 MPa [range, 6-18 MPa; 95% CI, 8-12 MPa] versus 8 MPa [range, 6-13 MPa; 95% CI, 7-9 MPa], p = 0.137). Deviation from acetabular sphericity positively correlated with preoperative peak contact stress (R = 0.326, p = 0.002) and was greater in the older cohort (0.9 mm [range, 0.8-1.5 mm; 95% CI, 0.8-1.0 mm] versus 0.8 mm [range, 0.6-0.9 mm; 95% CI, 0.7-0.9 mm], p = 0.014). Peak preoperative contact stress did not correlate with preoperative VAS pain score (R = 0.072, p = 0.229), and no correlation was found between change in peak contact stress and change in VAS score (R = 0.019, p = 0.280). CONCLUSIONS: Patients over the age of 30 years with dysplasia had less spherical acetabula and higher predicted preoperative contact stress than those younger than 30 years of age. Future studies with larger numbers of patients and longer term functional outcomes will be needed to determine the role of altered mechanics in the long-term success of PAO varying with patient age. CLINICAL RELEVANCE: These findings suggest that long-term exposure to abnormal joint loading may have deleterious effects on the hip geometry and may render the joint less amenable to joint preservation procedures. Given the lack of a direct relationship between mechanics and pain, orthopaedic surgeons should be particularly critical when evaluating three-dimensional dysplastic hip morphology in patients older than 30 years of age to ensure beneficial joint reorientation.

Unaddressed Cam Deformity Is Associated with Elevated Joint Contact Stress After Periacetabular Osteotomy.

BACKGROUND: Femoral cam deformity is frequently present in patients with acetabular dysplasia. Computational modeling can be used to identify how this deformity affects joint mechanics. Our purpose was to identify the relationship between cam deformity and joint contact stress after periacetabular osteotomy (PAO). We hypothesized that cam deformity is associated with an increase in peak joint contact stress after PAO. METHODS: This was a retrospective review of patients treated for hip dysplasia with PAO without femoral osteochondroplasty. Patient-specific hip models created from preoperative and postoperative computed tomography (CT) scans were evaluated using discrete element analysis to determine maximum joint contact stress after PAO. Twenty hips with a postoperative increase in maximum contact stress were compared with 20 that demonstrated decreased maximum contact stress. Hips were assessed for cam deformity on cross-sectional imaging. Radiographic measures of acetabular dysplasia before and after PAO were assessed and compared with the change in maximum contact stress after PAO. RESULTS: There was a moderate relationship between the change in maximum contact stress and the α angle (r = 0.31; p = 0.04), and the average α angle in the hips with increased maximum contact stress was significantly different from that in the hips with decreased joint contact stress (51° ± 11.4° versus 42° ± 5.1°; p = 0.04). All 6 hips with an α angle of >60° demonstrated increased joint contact stress. CONCLUSIONS: Cam deformity is common in patients with hip dysplasia. In our study, α angles of >60° were associated with increased postoperative joint contact stress. The α angle should be assessed preoperatively, and deformity should be addressed for optimal joint mechanics after PAO. CLINICAL RELEVANCE: A reduction in joint contact stress is a proposed mechanism for the increased joint longevity following periacetabular osteotomy for hip dysplasia. Impingement from abnormal femoral offset negatively impacts clinical outcome, but this finding has not been evaluated from a biomechanical perspective previously and a threshold for performing femoral osteochondroplasty has not been established previously. This study provides biomechanical evidence supporting surgical management of femoral cam deformity for an α angle of >60°.

Joint contact stresses calculated for acetabular dysplasia patients using discrete element analysis are significantly influenced by the applied gait pattern.

Gait modifications in acetabular dysplasia patients may influence cartilage contact stress patterns within the hip joint, with serious implications for clinical outcomes and the risk of developing osteoarthritis. The objective of this study was to understand how the gait pattern used to load computational models of dysplastic hips influences computed joint mechanics. Three-dimensional pre- and post-operative hip models of thirty patients previously treated for hip dysplasia with periacetabular osteotomy (PAO) were developed for performing discrete element analysis (DEA). Using DEA, contact stress patterns were calculated for each pre- and post-operative hip model when loaded with an instrumented total hip, a dysplastic, a matched control, and a normal gait pattern. DEA models loaded with the dysplastic and matched control gait patterns had significantly higher (p = 0.012 and p < 0.001) average pre-operative maximum contact stress than models loaded with the normal gait. Models loaded with the dysplastic and matched control gait patterns had nearly significantly higher (p = 0.051) and significantly higher (p = 0.008) average pre-operative contact stress, respectively, than models loaded with the instrumented hip gait. Following PAO, the average maximum contact stress for DEA models loaded with the dysplastic and matched control patterns decreased, which was significantly different (p < 0.001) from observed increases in maximum contact stress calculated when utilizing the instrumented hip and normal gait patterns. The correlation between change in DEA-computed maximum contact stress and the change in radiographic measurements of lateral center-edge angle were greatest (R2 = 0.330) when utilizing the dysplastic gait pattern. These results indicate that utilizing a dysplastic gait pattern to load DEA models may be a crucial element to capturing contact stress patterns most representative of this patient population.

Discrete element analysis is a valid method for computing joint contact stress in the hip before and after acetabular fracture.

Evaluation of abnormalities in joint contact stress that develop after inaccurate reduction of an acetabular fracture may provide a potential means for predicting the risk of developing post-traumatic osteoarthritis. Discrete element analysis (DEA) is a computational technique for calculating intra-articular contact stress distributions in a fraction of the time required to obtain the same information using the more commonly employed finite element analysis technique. The goal of this work was to validate the accuracy of DEA-computed contact stress against physical measurements of contact stress made in cadaveric hips using Tekscan sensors. Four static loading tests in a variety of poses from heel-strike to toe-off were performed in two different cadaveric hip specimens with the acetabulum intact and again with an intentionally malreduced posterior wall acetabular fracture. DEA-computed contact stress was compared on a point-by-point basis to stress measured from the physical experiments. There was good agreement between computed and measured contact stress over the entire contact area (correlation coefficients ranged from 0.88 to 0.99). DEA-computed peak contact stress was within an average of 0.5 MPa (range 0.2-0.8 MPa) of the Tekscan peak stress for intact hips, and within an average of 0.6 MPa (range 0-1.6 MPa) for fractured cases. DEA-computed contact areas were within an average of 33% of the Tekscan-measured areas (range: 1.4-60%). These results indicate that the DEA methodology is a valid method for accurately estimating contact stress in both intact and fractured hips.